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C3 . Abstract 

We propose a novel derivation of the non-local heat kernel expansion, first studied by Barvin- 

sky, Vilkovisky and Avramidi, based on simple diagrammatic equations satisfied by the heat 

kernel. For Laplace-type differential operators we obtain the explicit form of the non-local heat 

kernel form factors to second order in the curvatures. Our method can be generalized easily to 

O ■ the derivation of the non-local heat kernel expansion of a wide class of differential operators. 

c^ ■ 1 Introduction 

The heat kernel plays a central role in many areas of mathematics and physics. At present 
date two main expansion schemes have been developed to calculate the heat kernel for general 
^N , backgrounds. The first is the early time or local expansion, which is also known as the Schwinger- 

deWitt technique [TJ [5J [3]. The second one is the non-local expansion studied by Barvinsky, 
Vilkovisky and Avramidi [HIS]. This latter case can be seen as an expansion in the curvatures. 
(^ ' In quantum field theory the local heat kernel expansion is usually employed to compute ultra- 

CN ' violet divergences and anomalies, since these are directly related to local heat kernel coefficients. 

rf-\ I The non-local expansion is used instead to compute, in a fully covariant framework, the finite 

(^ . part of the effective action [4j [5] . It also has many applications to the effective average action 

CN I [5] and to the spectral action [7]. 

The aim of the paper is to introduce a new diagrammatic technique for the computation of the 

heat kernel, based on the combination of a properly defined vertex-expansion and momentum- 

J^ . space rules. These momentum-space rules can be derived from a single functional that we call 

''T^ ' "Laplacian action" , that serves as a device for the computation of the heat kernel. We will show 

C^ , how this method can be used to give a new and independent derivation of the non-local heat 

kernel expansion first derived in [3]. We will consider the non-local expansion for second order 

covariant Laplacians on a general boundaryless d-dimensional manifold with arbitrary gauge 

connection. A similar method has been proposed in [5], but our presentation has the advantage 

of bringing all the strength of Feynman-diagrams' reduction techniques into the computation of 

the heat kernel expansion. 

For a review of the more mathematical and geometrical aspects of the heat kernel see [S], 
while for a physicist's perspective see [THl E]. A fully covariant scheme of the expansion of the 
heat kernel and some related operator may be found in |12j . 

2 Basic definitions 

Let us work on a general boundaryless d-dimensional manifold, equipped with a Riemannian 
metric g^i, and inverse g^'^ . We call the unique torsionless connection compatible with the metric 
F^ and the associated covariant derivative V^, so that V f^g^p = 0. We also introduce a general 
vector bundle with connection A^ and define the full covariant derivative to be D^ = d^ + A^ . 
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We do not specify further the structure of the bundle, it suffices to know that the connection A^ 
is expanded in terms of a basis of generators T" as ^^ ~ A'^T"' , where repeated indices always 
imply summation. Notice that the connection A^ may coincide with the Levi-Civita connection 
r^ in some situations, for example when the vector bundle is the tangent bundle itself. 
The curvature of the connection is given by 

= d^,A,-d,Af,^[A^,A,], (2.1) 

and has components F^i, = F" T". In the simple case in which ()2.ip acts on the tangent vector 
bundle, we have by definition that f2^i/ = TZfj,^, where 7?,^i/ is the Riemann curvature 2-form. In 
components, the Riemann curvature 2-forin is the Riemann tensor R^,^" p. The Ricci tensor is 
defined as i?^,y = Rafi^iy, while the scalar curvature is defined as i? = g^'^ R^,y. 

We define a general Laplace-type operator to be a second order differential operator of the 
form 

A = -.g^'^D^D^ + U = -1)2 + U , (2.2) 

where an cndomorphism U over the vector bundle has been introduced. From a physicist per- 
spective, we may call U the potential term of the Laplace operator and its components are 
defined as U = UT". For short, A will often be called simply Laplacian. 

The heat kernel K^{x,y) is a bi-tensor of density weight 1/2 defined through the following 
partial differential equation and boundary condition: 

{ds + A,.)K%x,y)^0 K''{x,v)^5{x-v). (2.3) 

Notice that, when the Laplacian acts on tensors of different weight, also the coordinate repre- 
sentation changes. Wc temporarily adopt two different symbols to stress this feature. For any 
given tensor -0 of weight w = 1/2, wc have that 

A.,^=i/2^(a;) = g^'\x)A^.^^=o{9-^/\x)^{x)) . (2.4) 

The operators A^. „,^]^ /2 and A^, .u,=o can be understood as the components of the same abstract 
operator A, when acting on tensors with different weight. An important property is that the 
two share the spectrum. We could have, in principle, adopted a different weight normalization 
for (|2.3p . but the one wc choose turns out to be particularly convenient in the following [lOj . 
The heat kernel equation (j2.3p can be interpreted as describing a continuous diffusion process 
on the manifold. In this picture, the diffusing particles are all concentrated at the origin x = y, 
when the heat kernel proper-time parameter s is zero. The heat kernel proper-time s is related 
to the real time t and the diffusion constant V through the relation s = Vt. 
Equation p.3p is formally solved as 



K%x,y) = e'^^-S{x-y), (2.5) 

so that its trace is formally the following: 

Tri^" = Tre-'"^==(5(x - y)]^^^ = Tre"""^ . (2.6) 

In order to give a precise meaning to the rather formal expressions (|2.5p and (|2.6p . we want to 
introduce a basis of orthonormal eigenfunctions for A. Special care has to be taken when dealing 
with a non-fiat background where the initial condition in p.3p requires a particular choice of the 
eigenfunction basis. Wc introduce a basis 0„ of eigenmodes of weight w = 1/2 for the Laplacian 

A^(l)n{x) = XnMx) , (2.7) 

that is normalized as 

d''-x(j)n{x)(j)jn{x) ^ 6mn ^ (?!)„ (a;)0„ (j/) = (5(x - y) . (2.8) 



Equation (j2.5p is written rigorously using the basis (/)„ as 

i^^(x,2/) = ^e-^^"0„(x)0„(y), (2.9) 



where the initial condition of (j2.3p is satisfied thanks to 

if°(x, y)^Y. M^)My) = S{x - y) . (2.10) 

n 

Tracing ^2.9\i and using the orthogonality relations in (|2.8p . we obtain 

TrK' =ii I d'^xK'ix.x) =^e-'^- /"d''x0„(x)(^„(a;) = ^ e""^" . (2.11) 

Note that in (|2.8p there is no factor ^ in the integrals because of its weight normalization. This 
will turn out to be particularly useful later on, because it decouples the normalization of the 
basis from the metric dependence of the heat kernel (|2.5p . 

As a final remark of the section, we want to stress one of the many useful applications of 
the heat kernel. It is used in theoretical physics to give a precise meaning to and to evaluate 
functional traces. In fact, every trace of an arbitrary positive function of the Laplacian /i(A) is 
related to the heat kernel by means of a Laplace-transform 

Tr/i(A) = / dsh{s)T:Ye-''^ , (2.12) 

Jo 

where h[s) is the inverse-Laplace transform of h{x). Therefore, in order to compute such a 
functional trace, one simply needs to know the expansion of the trace of the heat kernel to any 
desired accuracy. The local and non-local expansion schemes mentioned in section [T] represent 
the two main expansion schemes with which (j2.12p has been successfully applied to quantum 
field theory computations. 

3 Perturbative expansion of the heat kernel 

In this section we develop the perturbative expansion for the computation of the heat kernel, 
where the covariant Laplacian is decomposed as the sum of a non- interacting Laplacian —0^ and 
of an interaction V in the following way: 

A=-52 + y. (3.1) 

The potential V contains U, all terms proportional to the gauge connection A^ and all terms 
obtained by expanding the metric g^v = 5^i, + h^i, in powers of the fluctuation h^v 

We explicitly give two simple examples of the decomposition (|3.1[) . The flat space limit is 



A=-92 + u (3.2) 

where simply 1/ = U; the gauge Laplacian where V contains all terms that vanish for Af^, ~ 0, 
A = - (a^ + A^) (9'' + A^) = -d^ -2A^a^ - d^A^ - A^^A^ . (3.3) 



In both cases there is no distinction between different density weights. 

We derive now the perturbative expansion in Y of the heat kernel. First we need to compute 
the heat kernel Kq of the operator —9^, around which we will perform the perturbative 
expansion. Here and in the following we use the compact notation K^ = K''{x,y) and S^y = 



^''''(x — y). We also define J^ = J d'^x and J ^ J ,2-kY ■ Analogously to (|2.3p . wc require that 
it satisfies the following equation with boundary condition: 

{ds - dl) Kl^y - Kl,.y - S,y . (3.4) 



The differential equation (|3.4p is easily solved in momentum space. If we Fourier transform the 
heat kernel 

Kl.y = / K^.„' e-'(-«+^'') , (3.5) 

Jqq' 

then p.4p becomes simply 

{ds + q') Klgg, = <,,' = K+,' ■ (3-6) 

The solution of (|3.6p is trivially seen to be 

K,,' = '5,+,,e-^«' . (3.7) 

Transforming this solution back to coordinate space, completing the square and performing the 
Gaussian integral gives the following result: 



1 (x-y)^ 

(47rs)rf/2 



^0%. - TIZT^--^^ ■ (3-8) 



This is the fundamental solution around which we will construct the perturbative expansion. 
To begin with, we first note that the heat kernel p.8p satisfies the composition rule 

^0,xy — / ^0,xz^O,zy ■ (3.9) 

J z 

We define the operator [/^^ = J^ KQ^^Kly. Using (|2.3p and p.6p . we find it satisfies 

J Z 

J z 

^0,xz^z^zy 

^0,XZ^zKQ^^^Uyjy . (3.10) 

We use Dyson's series to solve equation p.lOp in a compact form 

c^x, = rexp |- £ dt ^ A'o:L^.^L,} , (3-11) 

where the exponential is time-ordered with respect to the parameter s. Substituting p.lip into 
(|3.10p and using (|3.9p . we immediately find that 

Ky = I K^zTexp [-J^dtJ^ A^L^-^L^,} ■ (3.12) 

Rescaling the integration variable in (|3.12l) as i ^- t/s and using l\'S.9\i . we obtain the final formula 
for the perturbative expansion of the heat kernel 



Ky = [k^,xzTcxp\-s f dt f K-llV^K^l^y 

Jz L JO Jw 



^O.xy ^ 



dt /<irv,ifo^*. 



W dtj dt, K^':l;''W.K[t"''Wu.K^%y + OiV'), (3.13) 



that represents the core of the covariant perturbation theory of the heat kernel [4] . The finite- 
ness of each term in this expansion is guaranteed if V is smooth and the covergence has been 
extensively studied in [4] . The terms arc finite even if V has a discontinuity of a derivative or a 
delta- function singularity [13]. As long as U, A^ are such that V has the above properties and 
g^i, is smooth, the integrals in (j3.13p are convergent. 

The expansion p.l3p is conveniently represented graphically 



K' = 



-^f- 



+ s^ 



^^-^ h 



(3.14) 



where a continuous line represents a Ko^^y factor, while with each cross is associated an insertion 
of the interaction potential 14 . The general term of the expansion p.l3p reads 



(-s)" / dti--- f" ' dt„ 



t-s(l-tl)T/ 7>'S(tl-t2) TA T^St„ 



(3.15) 



We now use the expansion p.l3p to derive the perturbative expansion for the trace of the heat 
kernel. To do this we simply trace (|3.13p . This gives, in graphical form, the following expansion 
for the heat kernel trace: 




Tr A'^ 



and the general term in this expansion is now of the form 




(3.16) 



(-s)"tr / dh 
Jo 



dtr. 



K, 



s(l-ti+t„) 



21 \J.ZiZ2 



■■K, 



s{tn-2 — ' 
0,Z„-2Zn 



v.. 



(3.17) 



Note that in p.l7p we used the cyclicity of the trace to combine last and first fiat heat kernels 

of (ixra . 

Let $ denote collectively the set of fields $ = {U, A^, /i^i,}. It should be clear that both 
the trace of the heat kernel p.l6p and the heat kernel p.l4p are functionals of $. We define a 
n-point function of the trace of the heat kernel p.l6p to be 



S"TtK' 



J$(xi)...<5$(x„)' 



(3.18) 



while a similar definition can be adopted for the heat kernel as well. 

To make profit of the perturbative expansions p.l4p and p.l6p we derived for the heat kernel 
and its trace, we need to devise a method to actually compute the various contributions. The 
strategy that we propose in this paper is simple but effective: derive equations for the n-point 
functions of the heat kernel by taking functional derivatives of the expansion p.l3p with respect 
to \J, A^jhfii, and then set them to zero. We then extract all possible information from each 
equation so derived, by considering their Fourier transform in momentum space. The equations 
for the n-point functions of the heat kernel trace read as follows: 



TiK' 






*=o 



*=0 





S^TrK'' 



(5$(a;)(5$(y) 




(3.19) 




and so on. These relations will be used, as we show in the next section, to compute the heat 
kernel expansion order by order in the curvatures. The term curvatures is used in a generalized 
sense. We include in the set of curvatures U, r2^,y and all tensors constructed with R^u"^ p. In 
particular we will use the equation for the two-point function to determine the non-local form 
factors at second order in the curvatures. 

Our aim now is to obtain the momentum space rules to compute the various diagrams present 
in the equations for rt-point functions of the heat kernel expansion (j3.19p . It is useful, first, to 
define a functional, that we will call "Laplacian action" and that serves as a simple device 
to compute the vertices appearing in the expansion as Feynman rules. We first introduce an 
auxiliary field 4> taking values in the internal space and having density weight w = Q. The 
Laplacian action L[(j), $] is defined to be a quadratic action in 0, whose Hessian is the Laplacian 
(|2.2p with the correct density weight. This notion carries over straightforwardly to the case 
of more general differential operators of non-Lapalcian type. The propagation of the field (j) is 
directly related to the propagation of the heat kernel degrees of freedom. By construction, the 
second derivative of L[0, $] with respect to (j) is 

L(2;0) [^. <i,]^,^ = 5(x)i/4 A[<i>],, g{yr^'' , (3.20) 

where A[$] is the Laplacian appearing in (|2.2p . The dependence of the Laplacian on $ was 
made explicit for a better understanding. In p.20p two powers of the determinant of the metric 
are present, so that A[$]2;j, are the components of the Laplacian when acting on tensors of 
density weight w = as we are correctly densitizing the formula to match with the basis (|2.8|) 
we adopted. Therefore, from p.20|) we obtain 

= lJd''x^[g^''D^{^g-^^^)DM9-'^^)+9-'^^V<l>']- (3-21) 

We expand (|3.20p in a functional Taylor series in $, 

4f ^[0;$] ^ 42-°)[0;0] +^4',^i)[0;0]$. + ^ J^^LifJjO;0]^^<i>^ + 0{^') (3.22) 

to observe that Lxy [0;0] = —d^Sxy is precisely the flat space Laplacian around which we are 
expanding. We can thus write the interaction potential defined in p.ip as 



VJ.y = f L^y'AO; 0]$. + i / 4',1[0; 0]$,$^ + 0(<i>3) . (3.23) 

J Z "^ J zw 

This shows that the vertices in the diagrams p.lQp are precisely those constructed from the 
Laplacian action ()3.21|) . The Laplacian action generates all the interactions needed for the 
perturbative expansion of the heat kernel and therefore has the role of the generating function 
of our method [13] ■ 

We can now state the Feynman rules for the computations of p.l9p in momentum space. 
Bolder lines represent external $ lines, that are associated with functional derivatives w.r.t. 
$ in (|3.19p . Lighter lines instead represent the propagation of the auxiliary field </>, that has 
the flat space heat kernel as propagator. The arrows point in the direction of the momentum 
(when giving a vertex we always adopt the convention that all momenta are incoming). The 
composition of propagators is derived from (J3T3]) for any number of insertions by taking the 



Fourier transform of p.l5|) and chopping the V insertions. It is represented as 
Pi P2 Pn 

I 1 = jyhj\t,...£'^' dt^_,e-<'-'^^r>l^-sit.-t.)pl.. 



1 rti rtn-2 

dti 
Jo 

g-^(t„-2-t„-i)pf,_ig-^t„-iP^ ^ (324) 



where bold hncs arc placed in correspondence of any chopped V insertion for convenience. The 
V insertions, instead, appear as vertices. These are obtained by taking functional derivatives of 
(|3.2ip . evaluating the result in momentum space. They are given by 



Pi P2 



P3 ^ I ^ Pn 

Vi .. . 



= L^''"^')[0,0]p,.p„...,p„ (3.25) 



An integration over closed loop momenta has to be performed. In the next section we show how 
to apply these diagrammatic rules to the calculation of the heat kernel trace. 

The existence of Feynman rules for the diagrams p.24p and p.25p requires the presence of 
asymptotic momentum states in the manifold under consideration. This, in turn, implies that 
the manifold has to be asymptotically M'^. In a more general situation, it is still possible to 
develop the construction of the expansion of the heat kernel p.l3p around an operator that is 
more general than —d^, once a full basis of eigenstates of the latter is known. This general 
treatment however would go beyond the scope of this work. 

4 Calculating the heat kernel trace 

In this section we give the computation of all the possible graphs appearing in p.l9p using the 
rules (|3.24p and (|3.25p . The results have to be compared with an appropriate ansatz for the 
curvature expansion of the heat kernel. The most general expansion of the trace of the heat 
kernel (|2.5p up to second order in the curvatures U, Sl^j^ and all tensors constructed from R^^"' p 

is a 

+\RfR{sU)R + i?/flc,(sn)U + U/y(sn)U + 0^„/o(sn)r!^''] +0(7^3)} . (4.1) 

where tr stands for the trace over the internal vector bundle and □ = —D^ is the Laplacian with- 
out the endomorphism term. We refer the reader to appendix |^ for some mathematical details 
used in the contruction of (|4.ip . The expansion (|4.ip has the name of covariant perturbation 
theory of the heat kernel and its validity obeys the condition 

VVTe > -R? , (4.2) 

where TZ is any element of the set {i?^ti/, i?, U, il^^} [3]. We complement the condition (|4.2I) by 
requiring that the manifold under consideration is asymptotically flat, in order for momentum 
space rules p.24p and p.25p to exist. 

We realize that there are three constants to be computed, go; ffc/,o and gRfi- Further, at 
quadratic order in the curvatures, five so-called form factors fRtdx), fR{x), fRu{x), fu{x) and 
fn{x) appear [J. These represent an important departure from the standard Seeley-deWitt 
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method to compute the heat kernel in local form (based on an expansion in s), although it 
is possible to understand them as the resummation of infinite number of terms of the latter 
[5]. Each order in the curvatures of (|4.ip has been normalized by an appropriate power of s, 
so that all the constants and all the form factors are dimensionless if a canonical dimension is 
given to the covariant derivatives. It follows naturally, then, that the functions have to have the 
dimensionless product sD as argument. 

Notice that, since our manifold of choice is boundaryless all terms of the form n"i? and n"U 
for n > 1 are not present in (|4.1|) . We will come back to these terms in appendix IB] 



4.1 Endomorphism contribution to the heat kernel trace 



We start computing all the terms in the expansion of p.l9[) involving solely the endomorphism 
U as external line. We include in this set of diagrams also the diagram with no external lines at 
all, for simplicity. We need all the vertices with up to two external U lines, that in components 
are 




k-t, a 



h 



5L 



bkMkJvt. 



-[0,0] = r'^ 



(4.3) 



fc.-!,a 



h,b 



Ofcit 



SL 



hjvtsvi 



■[0,0] =0. 



(4.4) 



We adopted dashed lines for the external endomorphism U. 

The first diagram we compute is the single loop with no external lines. According to p.l9p 
and using (j3.24p in the simplest case of n = 0, in momentum space it is 



TtK" 



$=0 




trl = 



(27r) 



d'^qe-'"^'tTl 



trl 



(47rs)'^/2 



(4.5) 



This has to be compared with the ansatz (|4.ip in the collective limit $ = and therefore uniquely 
determines the constant go = 1. 

Then we compute the endomorphism 1-point function. Due to momentum conservation, the 
incoming momentum of the endomorphism line must be zero. We obtain 



STtK" 



sv- 



$=0 




n)d J 



(27r) 

1 
*(47rs)''/' 



trT" . 



(4.6) 



Again, we can compare (|4.6p with ()4.ip and determine gu.o = ^1- 

Finally, we want to obtain the endomorphism 2-point function. We observe that the vertex 
with two external U lines in (|4.3p is zero, therefore the expansion consists solely of one diagram 



S^TrK" 



SV^SVLp 



2s 



2 P 



<S>=0 




= 2s 



q + p 
1 rti 

^ ' dti dt2 I e 

Jo Jq 



-s(l-ti+t2)q^-s{ti-t2}{q+p)' 



trT'^r^) , (4.7) 



where momentum conservation has been used. It is convenient to introduce the new variable 
£, — ti — t2 and rewrite the exponential of (|4.7p as 

il-ti+h)q^ + ih-t2){q+pf = {q + ^p)^ + a^ - Op^ ■ (4.8) 

We can thus shift the momentum variable inside the integral of (|4.7p and change the parametric 
integrals into a ^-integration, that gives an additional factor of 1/2. These manipulations give 



S^TrK" 



f-/ ^^ pU 1^ _p 



$=0 



{AttsY/- 



:f{sp')tTT<'-T'>^ , 



(4.9) 



where we defined the function 



/(^) 



d^e 



-m-i)x 



(4.10) 



which is the basic form factor. The relation (|4.9I) has to be compared to the second functional 
derivative with respect to U of the ansatz (|4.ip , that is 



S^TtK' 



SV1,61J% 



2s' 



$=0 



[AttsY/- 



;fu{sp^)tTT<^'T'^ . 



(4.11) 



Equating (|4.9p and (|4.11|) we obtain 



fuix) 



:/(^) 



(4.12) 



4.2 Connection contribution to the heat kernel trace 



In this section we want to exploit all the terms of the expansion of p.l9p involving external lines 
of the connection A^. We adopt a curly notation for these external lines. The vertices we need 
in this subsection are 




fa, u 




k-j, II ki,v 



SL 



9fcl( 



^k^^K;,^ 



[0,0] = (fci-fc2)^T% 



SL 



bkMkJAt^.JAl.^^^ 



(4.13) 



[0, 0] = 2(5^"^ trT(°r^) , (4.14) 



where the parentheses imply symmetrization T'T'') = l/2(r"T'' + T'T"). 

It is possible to compute the 1-point function of the expansion p.l9p involving an external 
Afj^ line and show that it is zero. This agrees with the fact that no term linear in ^^ is present 
in the ansatz. Since this computation is not particularly instructive, we shall not show it here. 

We therefore turn our attention directly to the 2-point function. It consists of two diagrams 



S^TtK' 



"^p;/i<J^-p;jy 



2s^ 



$=0 





(4.15) 



that we compute separately. By the use of (|4.13p . the first diagram of (|4.15p gives 




dh dt2 / {2q + pY{2q + p)''e-'*(l-*i+*^)«'-'*(*i-*=)(«+P)'trr(T'') 



We perform the same manipulations that we did for the computation of (|4.7p and subsequently 
reduce the tensor structure inside the integration thanks to rotational invariance. The result is 




' ^-^^^"'^ '-■^^'^'%^p^YrT(^T''Y (4.16) 



{i-Ksy/^ 



fjsp^)^,, ^ l-.f{sp')^,^,^ 



sp^ 



where in the second line we used the fact that the basic form factor (|4.10p has the property 



dS,e 



-€(i-«)2; 



(1-20^ = --[/(•^)-l] 



(4.17) 



The second diagram of (|4.15p is much easier to compute and gives 




= 25^"' I e 
'q 



Sq^ ^j.rp{a.rpb) 



.^tJ-t^lj-J^io-J^b) 



Summing up (|4.16p and (|4.18p with the correct overall factors, we determine (|4.15p as 



^^TrA"' 



"^p;fj,^^-p;u 



<S>=Q 



(47rs)d/2^ T gp2 



(4.18) 



(4.19) 



where we introduced the transverse projector PJ^'^ . This and the longitudinal projectors are 
functions of the incoming momentum p^ and defined as 



P^" = 6^'' - ^ 



p.. ^ P^P-" 



p^ p^ 

The hessian of the ansatz (|4.ip is also transverse and is easily computed as 



^^TrAT^ 



^^p-p."^~p-u 



<S>=0 



{4:7:3^/' 



;P^PfMisp')ti-T'^''T''^ 



(4.20) 



(4.21) 



A trivial comparison of (|4.19p and (|4.2ip gives 



Mix) 



2x 



[/(x)-l] 



(4.22) 
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4.3 Metric contribution to the heat kernel trace 

In this subsection we complete the computation of the form factors appearing in the ansatz (|4.ip , 
including all the diagrams of p.l9p involving at least one external graviton line. We need to 
compute the vertices of p.2ip involving solely graviton external lines 



k[^kl^ -^(ki + k2)H'"' (4.23) 




2q(^5'')("/) + JpH^"'5°'f^ , (4.24) 



p, (jiv) -p, (q/3) 



where in the second one we do not give the most general dependence on the incoming momenta 
for simplicity. The first vertex of (|4.23|) is derived from p.2ip , by considering all the terms linear 
in /i^jy when g^^ ^ (5^^ + h^,^] similarly, the second vertex comes from all the terms quadratic 
in /ipj,. In this section we will also need the endomorphism vertices given in ()4.3p . In principle, 
we would also need vertices mixing the endomorphism and graviton lines, however, due to the 
density normalization of the basis (|2.8p we adopt, the endomorphism is decoupled from the 
metric in the Laplacian action p.2ip and therefore no mixed vertex is present. Further, since 
the ansatz (|4.ip does not possess any form-factor involving the curvature ll^i, and any of the 
tensors R, R^i, or Rfj,uai5, there is no need to report the vertex involving one external metric and 
one external connection lines. It is easy to see that, in this case, the 2-point function of (|3.19p is 
identically zero. In fact, the computation of the 2-point function of the heat kernel with external 
A^ and h^^, lines can be used as a check of the consistency of the ansatz. In the case of the trace 
of the heat kernel the result is zero, that implies that no form factor relates rj^j, with R or i?^^. 
Since this computation would not determine any of the form factors of the ansatz (|4.1|) . we shall 
not include it. 

Let us proceed computing the 1-point graviton function of (j3.19p 



Sho-niy 




<s>=o ^ , 

fri 

2(4..)^/' ■ ^'-''^ 



This has to be compared with the functional derivative of the ansatz (|4.ip in the same limit. It 
is easy to see that the result is consistent with the determination ^o = 1 ■ However (|4.25p cannot 
be used to compute g_R,o- This can be seen both from the ansatz and from (|4.25p . On the ansatz 
side, one shows easily that the first order expansion in the metric fluctuation ft,^^ around a flat 
background of the term of (|4.1|) corresponding to pfl.o is a total derivative and therefore vanishes 
upon integration for our boundary conditions. On the side of (|4.25p . we observe that momentum 
conservation forces the incoming momentum of the /i^^ external line to be zero and gfjQ cannot 
be detected since it would appear with a quadratic power of the incoming momentum. The 
correct way to compute gu^ in our formalism is by considering an ansatz for the heat kernel, 
rather one for the trace of the heat kernel. This will be shown in |B] For the moment, let us 
assume the well known value from the literature gafi = 1/6, that we are going to need in the 
following. 
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We now move to the computation of the form factors of (|4.1|) involving at least one curvature 
tensor constructed with the metric. The 2-point graviton function is 



S^TtK" 



^'^p;fiiy^'^—p;al3 



= 2s' 



*=o 





(4.26) 



The evaluation of the diagrams goes as in the previous cases, we find the following result: 



S^TrK' 



^'^p;/iz/^'^— p;a/3 



1 



$=0 



(47rs)'^/2 
1 



-l + lf{sp'] 



B 



.fj.i'.af3 



d-i 



-sp 



+ -^(4(d + 1) + 4(rf - l)sp^ + {d- l),sV)/(sp') 
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p; 



./ii/,a/3 



1 



5/11^, q/3 



1 pf^iy.alS , yd__}_( pt^v.afi 



^_pHV,al3 _|_ 



In (|4.27p we have introduced the following set of six generalized tensor projectors [TB] : 



P2;fj,u,af:) 
PS:,fj.i',ai3 



1 (pT pT , pT pT\ _ 1 pT pT 



— — P"^ P^ 

d-1 ^'^ "/^ 



1 



.pT pL 



1 



We list here some of their properties: 



pL pT 

-* /iiy-* a/3 






P,^ = P, , for i = 2,1,5', CT 

(Ps. + P,s)' = Ps + P^, 



(4.28) 



(4.29) 



that are useful when computing the components of (|4.27p . The Hessian of ()4.ip is also computed 



J^Tr/v" 



^'^p;fii-'^'^—p;a.p 



1 



*=0 



(47rs)'^/2 
+ 



.90 .9-R,0 2 I 1 2 4^ / 2 



B 



fj,iy.al3 



{d - 3).go , (d - 2)gjifi ^ 
: i ;; sp 



+ ^sV/ffic(sp2) + 2(d- l)sV/fl(sp') 



pfiv,aP 



90_pfj,u,al3 _ 90_pf_iu,aP , \/rf - IffO / pfj.u,ap pp. 



Scr 



-^^r'')}-(4.30) 



Now we use the fact that go = 1 as previously computed, and that gufi = 1/6 as will be shown 
in appendix |B] Equating (|4.27p and (|4.30p . we can solve the system of the coefficients of the 
projectors in terms of the form factors Jr{x) and frndx). We find that 



/ffic(x) = I- + 1[/(.t)-1] 

32 8a; 48a; 8a;^ 



(4.31) 



We finally want to compute the last form factor of the ansatz (|4.ip . namely funix)- To this 
end we first calculate the 2-point function of the heat kernel with external U and h^^, lines, that 
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is simplified by the fact tliat tlie endoniorpliisni does not couple to the metric explicit in the 
Laplacian action p.2ip . There is only one diagram contributing, 



S^TtK' 



Shp.^JV- 



2s' 



$=0 




(47rs)''/2 



-lp>i:^-\{2s + sV)P^\f{sp') 



trT" . 



that has to be compared with the corresponding Hessian of the ansatz (|4.ip 



S'^TvK" 



6hp.f,^S\Jtp 



<s>=o 



(47rs)''/2 



2 ^ 



gufis 



(4.32) 



sVfRu{sp^))Pr ta-r"- (4.33) 



Equating (|4.32p and ()4.33p one easily confirms that gu.o = — 1 and that 

huix) = -i/(x)-^[/(x)-l]. 



(4.34) 



5 Results for the non-local expansion 

In this section we briefly collect all the results for the non-local heat kernel expansion, analyze 
their properties and compare them with the existing literature. 
The expansion of the trace of the heat-kernel reads as follows 



Tr/\" 



(4 



-1— j d'^x^tv^l - sU + .sl| + s^ [lR^,fR,,{sn)R^'-' + lRfn{sn)R 



-RfRuisa)U + Vfu{sD)U + n,,,fnisD)n 



fj.l' 



o(7^3)} 



(5.1) 



where 7?. represents the endomorphism or any of the curvatures and □ = —D^. For the heat 
kernel form factors in ((??T|) . we collect the results of (|4.12p . (|4.22p . (|4.3ip and (|4.34p . These have 
been found to be 



fmdx) = ^ + ^[/(^)-l] 



fnix) = 


- 32^^^) + 8.^(^)^48. - 


1 

8x2 


fRuix) = 


= -\m~iifi^)-n 




fuix) = 


- ^/(-) 




fn{x) 


= -J-[/(x)-i], 





[/(^) - 1] 



(5.2) 



and all depend on the basic heat kernel form factor f{x) that is defined in terms of a parameter 
integral 



f{x)= / d^e 
Jo 



^^C(i-?) 



(5.3) 



Our findings are in complete agreement with |3], but given in a different curvature basis as we 
will show in subsection 15.31 We observe that (|5.ip is universal, in the sense that it depends on 
the dimensionality d of the manifold only through the overall factor {4'ks)~'^^'^. This is not a 
feature that can be expected from the heat kernel of any operator, but rather is a property of 
considering a Laplace-like operator. A plot of the basic form factor (|5.3p is given in figure [TJ 
Using in (|5.2p the Taylor expansion of the basic form factor ([57 



/(.)^l-M^ + 0(.3), 



(5.4) 
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Figure 1: Plot of the basic form factor f{x). The large x behavior shows an exponential tail, that cannot 
be obtained with finitely many terms in a small- a; (Seeley-deWitt) expansion. 



jives the following short time expansion for the form factors 

^2 



fmcix 

Iru{x 
fu{x 
fn{x 



1 a; 

60 ^ 840 



15120 
lla;2 



+ 0{x^) 



120 336 30240 ^ ' 



1 X 

6 ^30 



280 

„2 



+ 0(x3) 



i_^ + ^ + 0(^3) 

2 12 120 ^ '^ 



1 a; 

12 ^ 120 



1680 



+ 0(x3). 



(5.5) 



Inserting (j5.5p into (|5.ip it is possible to recover the results of the local expansion of the heat- 
kernel, for example those obtained with the Seeley-deWitt method. 

However, having access to the full form factors (|5.2p . we can study the late time behavior of 
the heat kernel as well. This is of interest, since by definition it is not accessible from the local 
expansion, unless the resummation of infinitely many terms in performed. From the late time 
expansion of the basic form factor 



fix) = - + 



O 



(5.6) 



we find 



Iric{x) 
fnix) 

fnuix) 
fu{x) 



O 



1 1 

6.T x"^ 

1 1 f I 

O 



1 



.T" 



2 8 ^f I 

X'^ x^ 

12 ( \ 



X x^ 



x-^ 



14 



Mx) 



1 

2x 



1 
2^ 



O 



(5.7) 



For some application of the form factors (j5.2p and the corresponding expansions (j5.5p and 
(j5.7p . we refer to [H [5]. Some comment on very important apphcations of (|5.2p in quantum 
field theory computations will be given in subsection 15.21 when a different curvature basis will 
be adopted. 



5.1 The two dimensional limit 

The two dimensional limit of (|5.ip is rather interesting and deserves a separate study. For d = 2 
the relation Rf^^, = ^QfiuR holds and therefore the form factors fR{x) and frndx) combine. 
Neglecting endomorphism U and connection, we have that 



TrisT" 



d=2 



^ I d^x^til [l + s| + s^Rh.2d(.sn)R + OiR^) 



where there is now only one form factor fR2d{x) defined as 



lR2d{x) 



1 



I{x) 



1 



32'' ' ' 16x 
The late time expansion of (|5.9p is found to be 

2 
fR2d{x) = ^ + 



Wix) - 1] 



8^t^(^)-^] 



12 



o 



(5.8) 



(5.9) 



(5.10) 



we observe that in (|5.10p both the coefficients of x ^ and x ^ are zero [4] . 

5.2 The Weyl basis 

It is possible to change the basis {Rfj.u, R, U, i7^i,} of the curvature invariants in (|5.ip to the set 
{Cfiuap,R,^,^fiu}, where C^^ap is the Weyl tensor [T5]. The expansion (|5.ip becomes 



TvK' 



i'^x^gtrj] 



6 



+iRfRMsisn)R + RfRuisn)V + U/c/(sn)u 



lC^.c.pfc{sa)C' 



fii^ajS 



~n^,fnisa)Q^"'^ +0(7^•■^)} 



(5.11) 



We defined two new form factors 

fc{x) 

IrUs {x) 



d-2 

d 
4(d^ 



■fRtcix) 
fRic{x) + fR{x) 



(5.12) 



The new basis spoils the universality property, because it induces an explicit dependence on d 
into the form factors. We refer to El for a more detailed derivation of (|5.12p . The short time 
expansions of the form factors (|5.1ip are fc{x) = (d — 2)/360(d — 3) + 0{x) and fRMs{x) = 
(3d-2)/240(d-l) + O(a;). 

The Weyl basis presented here is particularly useful in quantum field theory applications. In 
fact, the Hessian of the square of the Weyl tensor, when evaluated on flat space, is proportional to 
the projector P2 and is then related to the inverse momentum-space propagator of the irreducible 
spin-2 component of the fluctuation /i,j,y. On a general background, the form factor fc{x) is thus 
directly related to the full quantum propagator of "graviton" excitations |6]. 
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5.3 The BV basis 



We finally present the result (|5.ip in the basis of curvatures in which it was originally obtained 
by Barvinsky and Vilkovisky [J. We therefore introduce a new basis, that we refer to as the 
BV basis, where the trace of the heat kernel is expanded in powers of the set of generalized 
curvatures {Rfj,u, R, P, ^^v} with P = — U + lR/6. In the new basis (|5.ip becomes 



TiK" 



1 



i'^.TV^trj 



1 + sP + s^ 



(4^s)d/2 



ii?M-/i(sn)i?"" + ii?/2(sn)i? 



o(7^3)}, 



(5.13) 



and a new set of form factors fi{x) for i ~ 1, . . . , 5 has been introduced. These new functions 
are related to those of (j5.2p by a simple linear transformation 



h{x 



fmdx) 

fnix) + -^fu{x) + -fuR{x) 



1 



3 

fu{x) 

Mx) 



fu{x) - funix) 



(5.14) 



We refer to appendix 1X1 for a complete derivation of (|5.14l) . Substituting (j5.14p into (|5.2I) . it is 
possible to check by direct comparison that our results do coincide with those of [3]. 



6 Resolvent method 

In this section we introduce a slightly different formalism to handle the heat kernel trace which 
is based on the resolvent method. The method is meant to be a momentum space extension of 
|17j . capable of capturing the non-local features of the form factors computed in the previous 
sections. It is based on the following relation 



ide 

"2^' 



1 



(6.1) 



where is a complex variable and C is a countour in the complex plane having in its interior all 
the poles of the inverse function {x — 9)^^ of the right hand side. From now on, all the formulas 
have to be suitably analytically continued in 6 as much as possible. The relation (|6.ip is very 
useful since it allows to express the exponential of a differential operator, i.e. the heat kernel, in 
terms of the inverse of the differential operator A — 9 



TtK' 



id9 
"2^' 



-se 



Tr- 



1 



A-9 



Tr 



1 



L(2;o) 



(6.2) 



where we abbreviated Jg = L id9e ^^ /2tt and used the fact that A acts on tensors of weight 1/2 
inside the trace to write A = L^'^-'^'> [(j); $]. The inverse (A — 9)^^ is usually called resolvent of A. 
It is important that in (|6.2p the contour C contains in its interior all the poles of the resolvent. 
The application of functional derivatives with respect to $ on the right hand side of (|6.2p is 
obtained by virtue of the relation 6M~^ = —M~^5MM~^, that is valid for any operator M. In 
this framework, we observe that the appearance of the time ordered product in the perturbative 
expansion of the heat kernel given in equation p.l2p is due to the fact that the operator A and 
its variation 6 A do not commute in general, i.e. [A,(5A] ^ 0. At the price of introducing the 
auxiliary complex integration over 9 we gain an easier way to obtain equations for the n-point 
functions of the heat kernel trace. One simply has 



STrK' 
S<P{x) 



tr 



1 



L(2;0) 



L 



(2;i) 



yz 



L(2;0) 



wy 
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^^^^^ ^^^ trf^J^] L^:}}[^J^] L(-) 



S^x')S<i>{x) JeJy.^u. \L(m - 0j^^ — \Li^-^o) 6)^^^ ™-' \U^--^) 



vy 



1 \ ,(2;2) / 1 



^^^t^(^i(2;O)_0J^/.™..'U(2;O) J " (S-^) 



One can now set <& = in ()6.3p and move to momentum space in order to obtain a diagrammatics 
for the n-point functions of the heat kernel trace of section U that is shghtly different from 
p.l9p . This diagrammatic shares the vertices of the one developed in the previous sections. The 
difference, instead, lies in the fact that propagators arc closer to the standard ones of quantum 
field theory and, for this application, have complex Euclidean mass square —9. 

For the sake of the application, we re-obtain the results of subsection l4.2l in this new formalism. 
For simplicity, we restrict our attention to the simple case in which the Laplacian is A = —D^ 
with Df^ = d^ -\- Af^ on a flat four dimensional manifold with Euclidean metric (5^^. We also 
choose Af^ to be an abelian connection, to avoid the unnecessary complications of an internal 
non-commuting fiber structure. Since the two-point function is known to be transverse, we can 
directly compute it by contracting it with the transverse projector and obtain only its transverse 
component. There are two terms in the expansion 



S^TrK 




oAp-f^oA-p-i^ 



(6.4) 



We explicitly included a marked point in the expansion, because the propagator in the resolvent 
method does not satisfy the composition rule that is proper for the ordered product (|3.9p . This 
notation should make evident that there are three propagators in the first graph and two in the 
second. 

The transverse part of the first diagram of (|6.4p in momentum space is found to be 

/ !^ r '-a^da f d^ I l§!e-"(^-«-^^^(-^))e-- , (6.5) 

Jc 27r Jo 2 Jo Jq a 

where we used a standard relation of perturbation theory to combine the propagators 

-J— = f d£ ^- ^-T a^da f d£e-"(«^+(i-«'^) (6 6) 

A^B i''^{^A+{l-OBr 2 70 "^"io'^^' ' ^ ' 

and we shifted the momentum integration in such a way that the exponential is a function of 
q^ only. To obtain the final result one has to perform four integrations; it proves convenient to 
perform the a integration first. This integration has a convergence radius that depends on 9 and 
therefore the result has to be analytically continued as anticipated. A singularity in 9 appears 
explicitly in the denominator 

id9 r' f 2q^a'e-^o 

c 27r io A[9'-^+P'e(l-0]'' 

Now we use the theorem of residues to perform the 9 integration, where the only pole is 0* = 
q2 ^ p2^^-^ _ ^Y Further, we integrate over the momentum q and the parameter ^ to obtain 

-i dU — 27riRes,. ^ ^ = —-^-—2f{sp') . (6.8) 

The tadpole diagram of (|6.4p has only two propagators to exponentiate and these share the 
same momentum g^, therefore no parametric ^ integration will be present. Again, we will take 
into account only its transverse part, since we know that the final result has to be transverse 
in the incoming momentum. All integrations, with the exception of the parametric one, are 
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performed in the same order as for (|6.7|) and ()6.8p . using the residue theorem at 6'* = q^ ■ We 
obtain 






^ g^ 



R-cse, ■ 

2 

(47rs)'^/2 



(6.9) 



Needless to say, it is possible to explicitly check that the longitudinal parts of the two diagrams 
of (|6.4p cancel each other. 

We finally combine the contributions from the two diagrams (|6.8p and (j6.9p and compare 
with the transverse part of the Hessian (|4.ip shown in (|4.2ip . The result is 

/o(a^) = -^[/(x)-l], (6.10) 

that agrees with (|4.22p . It is interesting to note that, within the resolvent method, the parametric 
^-integration comes from the exponentiation of the propagators of (|6.5[) using (j6.6p . This is 
another manifestation of formula (|6.2p , that relates the heat kernel with the standard propagator 
of quantum field theory. 

We are not going to reproduce all the form factors with the resolvent method, since it would 
not be particularly instructive. However, it is important to stress that this method proves 
powerful when dealing with operators with a more complex structure than a Laplace-like one. 
In particular in |17| the heat kernel of a non-minimal operator is computed in its local form. 
Our results show how it is possible, with the momentum space technique, to address the issue 
of computing directly the non-local form factors of a non-minimal operator. We plan to come 
back to this issue in a future publication. 

7 Conclusions 

Wc proposed an algorithm to compute the non-local heat kernel expansion of a general Laplace- 
like operator on a Riemannian manifold, with a general connection and endomorphism term. 
The method is based on a Feynman-like set of diagrammatic rules and therefore is suited to 
incorporate all the strengths of the reduction algorithms developed for quantum field theoretical 
perturbation theory. 

The key idea of the method is that the diagrammatic rules are entirely derived from a single 
object that we called Laplacian action, that serves as a device to construct the perturbation 
theory of the heat kernel, through a generating function method. The analogy of the heat kernel 
with a pcrturbativc quantum field theory extends further: the existence of the Laplacian action 
allows to interpret the heat kernel as the generating functional of a quadratic and therefore solv- 
able theory of a field of density weight w — 1/2. This field possesses a non-standard propagator 
and is coupled to endomorphism, connection and metric as sources. 

All the computations we did were performed in momentum space. We used the momentum 
space technique to compute the expansion of the heat kernel up to second order in a set of 
generalized curvatures. Third order form factors can be derived along the same lines, although 
the results are comprehensively much more complex. The method is particularly powerful since 
it can be generalized easily to any given differential operator, including differential operators 
that are non- minimal and not of second order. 

We introduced the resolvent method as an alternative technique, that may be useful when 
dealing with more complex situations, in particular for non-minimal operators. The resolvent 
method pushes the analogy between the heat kernel and quantum field theory further. In fact 
the perturbation theory that is developed using the resolvent shares the standard propagator 
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of a field theory. We believe our method will prove very powerful in taming the very complex 
structure of the bare form factors of the heat kernel expansion, since it can borrow easily all the 
sophisticated techniques that have been developed for the computation of perturbativc quantum 
field theory. 

It is important to stress that our momentum space technique is very versatile and it should 
be evident that, with some trivial generalization, it can incorporate the computation of heat 
traces involving non-trivial operator insertions easily. As an example, one may want to compute 
the non-local form factors of the heat kernel trace in presence of covariant derivative insertions 
[18] . Therefore, our method can provide a non-local generalization of the results of [19] . 
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A Useful relations 

On a d-dimensional manifold one can define the Euler integrand 

which is a total derivative in d = 4. The Euler integrand is not a total derivative for any d. 
However, using Bianchi identities, one can prove [S] that for n > 

/ d-^x [Rf^^^pD^'R'"'"'^ - 4i?^^n"i?^'' + i?n"i?] = 0(7^3) ; (A.l) 

it follows that at order TZ^ there are only two independent form factors that can appear in (|4.ip . 
since one can reduce their number by one by using (jA.ip in a Taylor expansion of these. In 
general, we are free to eliminate the Riemann form factor 

f d''xR^,o.f!fme^isD)R^''''^ = f d''x{AR^,fR^e„^{sD)R^''' - Rfme,nisa)R} + 0(7^=^) 

by absorbing fRiem{x) into the definitions of fmcix) and /fl(a;). This leads to the ansatz (|4.ip . 
In a very similar way one can change the form factor basis from the one considered in this 
paper {R^i,,R, U, fi^^} to the Weyl basis {C^^ap.R, U, Vl^^}, where C^^af} is the Weyl conformal 
tensor, the square of which is 

4 2 

1^2 p r)fj,va(3 p PM^ _i Z P^ 

C - R,..,R - j^R,.R + (rf_i)(rf_2)^ ■ 

In this case we use the following relation that is valid up to third order in the curvatures 



-R 



to define the form factors: 



Jbms (x) = ^ Jmc [x) + !r (.t) : 
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R 



which are those introduced in equation ()5.12p of subsection 15.21 

Finahy we write the form factors in the basis used in subsection 15.31 The authors of [3] 
defined U = — P + IR/Q and computed the order 7i? part of the trace of the heat kernel in the 
form 

Jd''xti-[iR^,h{sa)R'^'' + iRf2{sa)R + Pf3isa)R + p/4(.sn)p + n^.f^isa)^^,} = 

j d''xtr[lR^,h{sU)R^''' + li?[/2(.sn) + \h{sU) + ^hisU) 

+ u[-/3(sn) - i/4(sn)] R + \JU{sU)\j + o^,/5(sn)f]^,} . 

We can match the right hand side of this expression with (|4.ip . In this way we find the following 
relations between the two sets of form factors: 

Iur{x) = -h{x) - -/4(a;) 
fu{x) = fi{x) 

fR{x) = f,{x) + lMx) + ^Uix) 
6 36 

fmc{x) = fi{x) 
h{x) = f5{x); 
by inverting these we obtain (j5.14p . 

B Heat kernel in the coincidence limit 



As pointed out in subsection l4.3l it is not possible to compute the coefficient qr^, that appears in 
(|4.ip . by considering the momentum space diagrammatic we developed for the trace of the heat 
kernel. The reason is that the momentum conservation of the 1-point graviton function forces 
the incoming momentum to be zero and therefore kills the first variation of (|4.ip with respect 
to the metric. In other words, the first variation of (j4.ip . when evaluated in fiat space, is a total 
derivative and therefore cancels because of boundary conditions. In this appendix we provide 
an extension of p.l9p . that allows to compute the heat kernel and therefore has access to the 
coefficient g^.o- For simplicity, we limit ourselves to the first order in the curvatures. 

We begin with providing the most general ansatz for the heat kernel, like we did for the trace 
of the heat kernel in (|4.ip . It is 

K%x,x) = .^^L/2 Vg{go 1 + s guisD) U + sl gnisO) R + 0(7^^)} , (B.l) 

where two new non-local form factors gu{x) and gR{x) have been introduced. The heat kernel 
of (jB.ip has been considered in the limit of coinciding coordinates y = x, also called coincidence 
limit. We can Taylor expand these two functions and integrate in x to compare (jB.ip with 
(|4.ip . to discover that only the first term of the expansions of the non-local form factors does 
survive the boundary conditions. The comparison also shows that the consistency between the 
two ansatzs strictly requires gu(fi) = gu.o ^nd gu{0) = gn.o- Our task is then to compute gu{x) 
and gii{x), that will uncover these two additional non-local structures, as well as give a check 
for the computation oi gufi and a determination of g/jo- 

The diagrammatic expansion of the heat kernel is easily obtained considering p.l3p . that 
has been represented graphically in p.l4p . Functional derivatives with respect to the fields $ 
act on the cross insertions of p.l4p and are represented by external lines. We are interested in 
the curvatures to the first order, thus the relevant diagram in the expansion is the 1-point heat 
kernel diagram. It reads 

(B,2) 



5^{x) 
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and in particular we want to consider $(a;) = U(a;) and $(a;) ~ h^^nix). The case in which 
$(x) = ^/j(a;) can be shown to give a zero result at the coincidence limit. 

We first compute the diagram with an external endomorphism U line in momentum space, 
that gives 



<SUj, 



V 

1 



qJn 
It is trivial to manipulate the integral in the same fashion of those of section 51 to get 

Comparing our computation with the first variation of (jB.ip . one obtains 

gjj{x) = -fix), (B.5) 

that has the correct limit gu{0) = ^/(O) = —1- Interestingly, also this form factor can be seen 
as the resummation of infinitely many terms of the local expansion of the heat kernel [5]. In 
particular, it contains all the information about the terms of the form □"U, that are suppressed 
in our original ansatz by the boundary conditions for any n > 0. As a further check, it is easy 
to see that the first order in the expansion gives the term 

s^ 1 , 

-DU, 



(47rs)d/2 6 



that agrees with the results of the local expansion [T] . 

We now proceed by computing gR{x), for which the 1-point function at $(x) = h^v{x) is 
needed. We obtain 



^'^p-jfll' 



V 



(^ + iv)./(V)pr + ^^r], (B.6) 



(47rs)d/2 

that has to be compared with the same quantity derived from (jB.l[) . As expected, the longitudinal 
component of this result does not depend on the incoming momentum, since its contribution 
arises solely from the first (volume) term of (|B.1[) . The complete comparison allows to determine 
gR{x) as 

9r{^) - -^ + ^(.T + 2)/(x), (B.7) 

and, in particular, we find 

9Rfi = 5fl(0) = \ . (B.8) 

D 

As for the previous case, the form factor gnix) can be understood as the resummation of infinitely 
many terms of the local expansion. In this case it resums all the possible terms of the form 0"R. 
We can therefore check, as an example, the first term of the expansion that contributes to the 
local heat kernel as 

'' 'or, 



{Ansy/^ 30 
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that agrees with the resuhs of the local algorithms [T] . 

It is interesting to note that the two form factors computed in this appendix, gu{x) and 
gnix), share the same basic form factor (|5.3p with the system (|5.5p . This is a general feature of 
the expansion of the heat kernel. Namely, the heat kernel at n-th order in the curvatures will 
share the same basic form factors of the trace of the heat kernel at (n + l)-th order as shown in 



22 



References 

[1] D. V. Vassilevich, Phys. Kept. 388 (2003) 279 (hep-th/0306138| . 

[2] A. O. Barvinsky and G. A. Vilkovisky, Phys. Kept. 119 (1985) 1. 

[3] G. A. Vilkovisky, CERN-TH.6392/92. 

[4] A. O. Barvinsky and G. A. Vilkovisky, Nucl. Phys. B 282 (1987) 163; A. O. Barvinsky 
and G. A. Vilkovisky, Nucl. Phys. B 333 (1990) 471; A. O. Barvinsky and G. A. Vilko- 
visky, Nucl. Phys. B 333, 512 (1990); A. O. Barvinsky, Y. V. Gusev, V. V. Zhytnikov and 
G. A. Vilkoviskv. [arXiv:0911. 11681 [hep-th] . 

[5] I. G. Avramidi, Phys. Lett. B 236 (1990) 443; I. G. Avramidi, Phys. Lett. B 238 (1990) 
92; I. G. Avramidi, Lect. Notes Phys. M 64, 1 (2000). 



[6] A. Codello, Annals Phys. 325, 1727 (2010), arXiv:1004.2171' [hep-th]. A. Satz, A. Codello 
and F. D. Mazzitelh, Phys. Rev. D 82, 084011 (2010), arXiv:1006.3808 [hep-th]; A. Codello, 
New J. Phys. 14, 015009 (2012^. larXw:1108. 19081 [gr-ac]. 

[7] B. lochum, C. Levy and D. Vassilevich. larXiv:1108. 37491 [hep-th]: B. lochum, C. Levy and 
D. Vassilevich, arXiv:120 1.6637. [math-ph]. 



I. G. Moss and W. Naylor, Class. Quant. Grav. 16, 2611 (1999), jarXiv:gr-qc/0101125 



[9] S. Rosenberg, Cambridge, UK: Cambridge Univ. Pr. (1997) 174 p. 
[10] V. Mukhanov and S. Winitzki, Cambridge, UK: Cambridge Univ. Pr. (2007) 273 p. 



[11] G. Esposito, Cambridge Lect. Notes Phys. 12, 1 (1998), hep-th/9704016 



[12] L. L. Salcedo, Eur. Phys. J. C 49, 831 (2007) |arXi^:hep-th/0606071 

[13] M. Bordag and D. V. Vassilevich, Phys. Rev. D 70 (2004) 045003 [hep-th/0404069] . 

[14] A. N. Vasiliev, Amsterdam, Netherlands: Gordon and Breach (1998) 312 p 



[15] A. O. Barvinsky, A. G. Mirzabekian and V. V. Zhytnikov, |gr-qc/9510037! A. G. Mirz 



abekian, G. A Vilkovisky and V. V. Zhytnikov, Phys. Lett. B 369 (1996) 215 
|hep-th/9510205| . 

[16] K. S. Stelle, Phys. Rev. D 16 (1977) 953; L. Griguolo and R. Percacci, Phys. Rev. D 52, 
5787 (1995), arX iv:hep-th/9504092, , 

[17] V. P. Gusynin and V. V. Kornyak, Nucl. lustrum. Meth. A 389 (1997) 365. 

[18] A. O. Barvinsky and Yu. V. Gusev, Sov. Phys. J. 34, 858 (1991); A. O. Barvinsky 
and Y. .V. Gusev, Class. Quant. Grav. 9, 383 (1992); A. Barvinsky and Y. Gusev, 



arXiv:gr-qc/9507026, Y. V. Gusev, Nucl. Phys. B 807, 566 (2009) a rXiv:081 1.1063, [hep- 



th]]. 
[19] Y. Decanini and A. Folacci, Phys. Rev. D 73 (2006) 044027, |arXiv:gr-qc/0511115| 



D. Anselmi and A. Benini, JHEP 0710 (2007) 099, arXiv:0704.2840' [hep-th]; D. Benedetti, 
K. Groh, P. F. Machado and F. Saueressig. .arXiv:1012.3081. [hep-th]: K. Groh, F. Saueressig 
and O. Zanusso. llarXiv:1112.4856l [math-ph]. 



23 



